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Disclosure of Invention 



"Significant Dust Enhancement Over Land and Ocean" 

Steven D. Miller 
Naval Research Laboratory 
Monterey, CA 93943-5502 

a) Purpose 

The purpose of the Land/Ocean Significant Dust Enhancement Product is to provide ah 
indication of the horizontal distribution of airborne dust over both land water surfaces for 
daytime scenes. The regions of dust are enhanced as shades of pink, depending, on 
opacity, against darker non-dust backgrounds. 

b) Background 

A recent study (R. Cantu; 'The Role of Weather in Major Naval Aviation Mishaps", MS 
Thesis, 106 pp., NPS Monterey) has indicated that poor visibility has factored' into 
approximately 54% of Navy aircraft "Class-A" mishaps, for an estimated annual loss of 
$51 million. Of these mishaps it was reported that 56% were preventable had better 
forecasting and/or observational tools been available. The current Dust Enhancement 
Product represents a paradigm shift in the way dust is observed from the space platform, 
and is directly applicable to mitigating this multi-million dollar liability. 

c) Description and Operation 

As presented, the Dust Enhancement Product is actually a combination of two 
independent enhancement algorithms — one applicable to water backgrounds and the 
other to land backgrounds. Miller and Lee (AMS 11* Satellite Meteorology Conference; 
October 2001, Madison WI, pp. 43-46) and Miller (NRL Formal Report; 2002, Accepted) 
provide a comprehensive discussion of the ocean-component of the Dust Enhancement 
Product as applied to the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) instrument 
(these publications have been attached to this application). It is based physically on the 
premise that silicate aerosols preferentially absorb blue light compared to green and red 
light, whereas cloud particles (liquid or ice) are more or less non-absorbing across this 
region and reflect all three components strongly. Therefore, a spectral difference 
between red and blue light will be much larger for dust than for clouds. In this way, data 
from a radiometer with spectral resolution capable of measuring the red, green, and blue 
components of reflected solar energy can be combined to further enhance and decouple 
the dust signal from clouds. The spectral difference is normalized by the magnitude of 
the reflectance such that additional sensitivity to thin (low reflectance) dust is achieved. 
The approach requires . correction of the multi-channel data to remove scattering 
contamination from the molecular atmosphere prior to computing the enhancement. 
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No literature discussing the over-land component of the current Dust Enhancement 
Product has been drafted to date, although a submission to the peer-review literature on 
the topic is planned. A PowerPoint presentation, providing an overview of the method as 
applied to the Moderate Resolution Imaging Spectroradiometer (MODIS, aboard the 
Earth Observing System Terra satellite) and examples, is available upon request. 
Because the land background possesses similar spectral properties to the airborne dust at 
visible/shortwave-infrared wavelengths, the simple ocean-component described above is 
insufficient for identifying dust over land. We must instead appeal to thermal contrast 
and spectral property differences, which require channels in the infrared part of the 
spectrum. With its 36 narrowband channels spanning the optical portion of the 
electromagnetic spectrum fro 0.4 to 14.4 micrometers, MODIS is well equipped to 
exploring these techniques from a single platform (although the Dust Enhancement 
Product may be achieved by combining registered data from a collection of independent 
platforms). 

Three factors contribute to the over-land Dust Enhancement Product, and are outlined as 
follows: During the day, cooler emissions from an elevated dust layer contrast against 
the warm emissions from a heated surface. The same temperature signature holds true 
for clouds, but the same spectral decoupling applied to the over-ocean enhancement 
applies also to clouds over land. Combining the temperature and shortwave-differences 
provides a means to detecting the optically thick portion of the dust plume. Optically thin 
regions of dust will not produce a strong infrared temperature depression. However, 
spectral differences in thermal (11.0-12.0 micrometers) transmission for dust are very 
useful in determining areas of thin dust. The signature is opposite in sign to that of thin 
cirrus. Combining these multispectral features together yields a tractable algorithm for 
the enhancement of dust over land. Although the land and water enhancement algorithms 
differ significantly, the land enhancement has been formulated such that dust over land is 
the same hue (shades of pink) as dust over water with minimal discontinuity in saturation 
across coastal boundaries. 

d) Advantages and New Features 

The primary advantage of the Dust Enhancement Product is its ability to call immediate 
attention to areas where there airborne dust resides, often in regions where even the true 
color imagery experiences difficulty in discerning. The enhancement highlights 
mesoscale circulations in optically thin dust, and is also useful for identifying volcanic 
plumes (e.g., Mt. Etna 7/26/01 eruption — image attached). The latter capability is of 
particular relevance to both military and civilian aviation. The algorithm is applicable to 
a wide range of ocean color instruments in the ocean-only mode (i.e., true color over land 
surfaces) or in the ocean/land mode provided infrared information is available either on 
the same platform or co-registered with temporally close imagery from an independent 
platform. 

e) Alternatives 

Alternative methods exist for enhancing dust over land and water. Panchromatic 
(black/white) imagery with scaling thresholds chosen to enhance a small dynamic range 
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of temperature or reflectance values has been applied to visible and infrared geostationary 
imagery. The resultant enhancement appearance, usually only of marginal quality, is 
highly dependent on the thresholds chosen and hence not suitable for operations. 
Combinations of single-channel visible and infrared data is a marked improvement over 
the scaling described above, but still suffers cloud/dust ambiguity for lack of multi- 
channel visible data (i.e., the basis for the over-ocean Dust Enhancement Product which 
requires blue and short wave infrared channels). The Empirical Orthogonal Function 
(EOF, or sometimes referred to as "principle component") method has been shown to do 
a good job of enhancing dust over land and water. However, results produced by this 
method are highly dependent on the variability of the scene and hence will provide an 
inconsistent enhancement depending on background and dimension of the region. 
Contrast reduction, which takes advantage of the blurring affect of dust over otherwise 
sharp terrestrial features, has been applied successfully for dust detection over land, but 
naturally fails over ocean and also struggles in regions where terrain is less variable (e.g.-, 
desert plains, where dust incidentally is most common). 

None of the above-mentioned methods enlists multi-channel visible data to take 
advantage of the inherent spectral variability of dust at these wavelengths. The unified 
land/ocean Dust Enhancement Product is a unique solution to the daytime dust detection 
problem. A comparison between the visible/infrared combination technique (applied 
here to Meteosat 5 data) and the current algorithm (applied to Terra-MODIS data) for a 
space and time co-located dust storm event in Southwest Asia is attached. The superior 
detection of dust by the current method over both land and water as demonstrated by 
these data are thought to be of immediate relevance to Navy METOC operations in 
several regions of the world that experience significant dust outbreaks (including 
Southwest Asia, Africa, the Mediterranean, and the West Pacific). 
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EXECUTIVE SUMMARY 

INTRODUCTION 

A new processing capability for true color, vegetation enhancement, and desert dust 
enhancement using imager data from the Sea-viewing Wide Field of View Sensor (SeaWiFS) has 
been developed for Naval meteorology/oceanography (MetOc) operations support. The data are 
captured via direct broadcast high-resolution picture transmission (HRPT) at Navy Regional 
Centers in Rota, Bahrain, and Yokosuka, and processed at the Naval Research Laboratory in 
Monterey. The raw data are calibrated, corrected for missing lines and clutter, corrected for 
molecular scatter contamination, and enhanced through multispectral combination to yield value 
added products. The processing has been automated completely such that products, generated 
upon receipt of data, are hosted upon a password protected website typically 60 to 90 minutes 
from time of initial capture. This report summarizes the SeaWiFS instrument capabilities, the 
protocol followed for automated near real-time processing, a physical basis for the NRL 
enhancements, and specific examples of the products. It closes with a glimpse of the potential 
utility of these products from the perspective of the warfighter. 

FINDINGS AND RECOMMENDATAIONS 

It has been demonstrated that the SeaWiFS instrument has significant potential to benefit 
the warfighter through its ability to provide high resolution (1.1 km) visible imagery of 
' atmospheric dust outbreaks with enhancements for optically thin regimes, as well as other 
meteorological/land-surface features of interest (e.g., clouds, fires, snow cover). Based on these 
findings, NRL and Fleet Numerical Meteorology and Oceanography Center (FNMOC) have 
collaborated to host a subset of these products upon the FNMOC Secure Internet (SIPRNET) 
beta-webpage for the purpose of i) making this new technology available to the warfighter, and ii) 
receiving feedback useful for improving its utility. It is recommended that the Navy continue to 
pursue the procurement of all telemetries having true color imaging capability to address the poor 
temporal resolution inherent to the polar orbiter platform. In this capacity, imagery from 
MODIS, OCM, and FY-1C are obvious candidates, with the former two requiring X-band 
receiving stations. 



3 



SeaWiFS True Color, Vegetation, and Dust Enhancement Processing at NRL 

Monterey 

INTRODUCTION 

The Satellite Meteorological Applications Section (Code 7541) at the Nayal Research 
Laboratory (NRL) in Monterey has developed true color, vegetation enhancement, and airborne 
sand/dust enhancement processing capabilities for the Sea-viewing Wide Field of View Sensor 
(SeaWiFS) instrument. These products are made in near real time (typically 60 to 90 minute 
latency from time of ingest at the receiving stations) using telemetries captured at Navy Regional 
Centers (NRCs) located in Spain (Rota), Bahrain, and Japan (Yokosuka), and populate a 
password-protected website for several predefined sectors within the general coverage regions of 
these centers. Presented here is an overview of the SeaWiFS sensor and its 
capabilities/limitations, a summary of the procedure, and the current processing architecture. 
Examples drawn from the archived database (automated processing of these data commenced 
August 8, 2001) are presented to illustrate their quality and potential operational utility. 

SeaWiFS— A Technical Overview 




Fig- 1 The SeaWiFS sensor (left) and SeaStar platform (right) 
Platform/Instrument Specifications 

The SeaWiFS instrument (Firestone and Hooker, 1996), launched aboard an extended 
Pegasus rocket on August 1, 1997, is an 8-channel whiskbroom scanning radiometer featuring 8 
narrowband (-20 nanometers (nm) wide) spectral channels centered at 412, 443, 490, 510, 555, 
670, 765, and 865 nm wavelengths. Carried aboard the SeaStar spacecraft (see Figure 1), 
SeaWiFS follows a sun-synchronous (98.3° inclination, 705 kilometer (km) altitude above mean 
sea level (AMSL)) orbit with a local noon descending equatorial crossing time. It has a local- 
area-coverage (LAC) swath width of 2801 km (corresponding to a maximum scan angle of 
58.3°) and a nominal spatial resolution of 1.1 km at nadir. Global Area Coverage (GAC) is also 
available at 4 km degraded resolution and a reduced 1502 km swath (data limited to 45° scan 
angle). The 10-bit digitized data are transmitted down to Earth at an L-band frequency of 
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1.70256 Gigahertz (GHz). The instrument scans across 1285 picture elements at a rate of 6 Hz. 
Typical High Resolution Picture Transmission (HRPT) station telemetries include typically 4000 
to 5000 scan lines of data (depending on local topographical/anthropogenic obstructions), 
translating to 11-13 minutes of capture time from the nominal initial-capture time stamp on the 
filename. 

The primary function of SeaWiFS is to detect and monitor ocean color. Specifically, it 
has been designed to retrieve the small component (-5% of total signal) of water leaving reflected 
radiance whose spectral signature provides useful information about ocean phytoplankton content 
(expressed in terms of chlorophyll-A retrievals). A varying along-track sensor tilt angle of +/- 
19.8° is used to avoid regions of sun glint (specular reflection of the solar disk upon the ocean 
surface) where chlorophyll-A retrievals cannot be performed. The information obtained from 
such retrievals is valuable to marine research and industry interests. Applications include 
guidance to the fishing industry, diver visibility information, surface currents, physical 
oceanography, general marine biology, and climate impacts on the marine ecosystem (e.g:, 
McClain et a/., [1995], and Gordon et aU [2001]). While atmospheric constituents such as cloud, 
dust, and aerosol are considered as noise to ocean color algorithms, the current research regards 
these components as the signal of interest. This report describes how SeaWiFS data are applied 
to enhance atmospheric and surface features/phenomena. 

Instrument Advantages 

The main advantage provided by SeaWiFS over optical spectrum instruments flown on 
the geostationary (e.g., the Geostationary Operational Environmental Satellite (GOES) 
constellation, Meteosat 5,7, Geostationary Meteorological Satellite (GMS), and the Indian 
National Satellite System (INSAT)) and other polar orbiter (e.g., National Oceanographic and 
Atmospheric Administration (NOAA), Defense Meteorological Satellite Program (DMSP), and 
many others) platforms for the current application is its multiple narrowband visible channels that 
provide a true color imaging capability. The Moderate-resolution Imaging Spectroradiometer 
(MODIS; King et al, [1992]) carried aboard Earth Observing Satellites (EOS) Terra and soon to 
be launched Aqua platforms also provides these channels and additional channels reaching into 
the infrared (36 channels between 0.4 and 14.4 Jim). At the time of this writing, MODIS data are 
only available to NRL over the areas of interest to this study on a research basis at 2-week latency 
from the Goddard Space Flight Center (GSFC) archives. MODIS telemetry is transmitted at X- 
band (7^8.5 GHz) whereas Navy facilities currently are limited to receipt of L-band (1.7-2.5 
GHz). Direct broadcast MODIS data have been made available via the internet from the Space 
Science and Engineering Center at the University of Madison, Wisconsin, at the following 
address: http://eosdb.ssec.wisc.edu/modisdirect/ . courtesy of Dr. Liam Gumley of the Cooperative 
Institute for Meteorological Satellite Studies (CIMSS). 

Instrument Limitations and Other Caveats 

The main deficiencies of SeaWiFS in terms of 24-hour operational utility are its lack of 
infrared channels, its limited temporal resolution, and proprietary nature. Infrared data are useful 
for nighttime imaging, cloud particle size estimation, and airborne dust enhancement over land 
surfaces, among other uses. At mid to low latitudes, only one to three passes are available near 
local noon — providing at best a "snapshot" of the synoptic-scale situation on any given day. 
Captures are increasingly frequent at the higher latitudes owing to the orbital geometry. Typical 
captures for the regional centers currently processed are on the order of 2 to 3 per day. Thus, 
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SeaWiFS products are most useful when juxtaposed with geostationary data that offer the 
temporal information required to infer details of the atmospheric circulation. 

Because SeaWiFS was bui lt> launched, and operated by Orbital Sciences Corporation 
(OSC), the data are formally available to the general public on a research-only basis and are 
provided at cost. The National Aeronautics and Space Administration (NASA) has contracted 
with OSC for the research use of a 5-year data set commencing in September 1997. Researchers 
who have registered as authorized SeaWiFS user (registration at 
http://seawifs.esfc.nasa.go v/SEAWIFS.htmn are granted permission to obtain the data from the 
NASA SeaWiFS archive, approved NASA Space Act Agreement ground station, or authorized 
individual. As authorized receiving stations, the NRCs providing near-real-time LAC data for 
this research are covered under the existing agreement, with the only criterion being that NRL 
does not disseminate the LAC-derived products to the general public in near real-time. 

Data Processing 

Several options exist for processing the SeaWiFS telemetry, depending on the format of 
the data. The formats currently received at NRL are NASA Level-0 (LO; 10-bit Data Capture 
Facility frame-formatter files), with Terascan Data Format (TDF) (pre-calibrated and navigated) 
upon special request. Data from the NRCs are in NASA LO format, and require several 
preprocessing steps to calibrate and navigate the data as described below. For reasons described 
further along, the current approach to creating SeaWiFS products is based on NASA LO data and 
is independent of the Terascan software package. Several approaches have been iterated upon in 
arriving at the current procedure including various blended combinations of software packages — 
some more appropriate for rapid deployment to remote operational centers and others geared 
toward in-house processing with all the conveniences of locally available third party software. 
While the purpose of this report is to summarize the current approach, alternative options and 
their associated data/software requirements will be touched upon in passing. 

Worth noting is that the Perl programming language (Wall et aL, [2000]), running on the 
Linux platform, was selected as the scripting language for automated processing of SeaWiFS 
data. Perl draws from the strengths of many popular scientific programming languages such as 
C++ and Fortran while enhancing on the capabilities of other standard shell scripting languages. 
The Perl scripts written in support of the current SeaWiFS processing perform the requisite 
secretarial operations of data retrieval, file and directory manipulation, execution and status 
monitoring of telemetry updating and data processing codes, dissemination of imagery products 
to the website, and general housekeeping. Perl remains the lone constant among all processing 
options discussed in this report. 

SeaDAS Processing Option 

One of the available SeaWiFS processing options is the GSFC SeaWiFS Data Analysis 
System (SeaDAS) software package (Fu et a/., [1998]), which was designed specifically for 
viewing and basic manipulation of SeaWiFS data. It uses an embedded runtime Interactive Data 
Language (IDL) license as its operating paradigm. The SeaDAS software is designed to operate 
on NASA LO, Level- 1 A (LI A; raw spacecraft and instrument telemetry retained with converted 
geo-location, instrument telemetry, and selected spacecraft telemetry appended), and Level- IB 
(LIB; sensor calibration applied to the data). LI A and LIB data are stored in Hierarchical Data 
Format (HDF). The software includes a graphical user interface enabling straightforward 
viewing and manipulation of the data, as well as command-line functions useful for incorporating 
into scripts as batch jobs. The SeaDAS package was used early on in this research to examine the 
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general utility of SeaWiFS data in terms of coverage and data quality before proceeding with 
customized application. While SeaDAS is not used in the current processing (in the interest of 
gaining a higher level of control over product display options), several functions were adopted 
from this package for conversion from NASA LO to LIB HDF files and data correction. SeaDAS 
remains an extremely usefiil stand-alone tool for SeaWiFS data visualization, and is available via 
the Internet at http://seadas.gsfc.nasa.gov . The website includes versions of SeaDAS for several 
operating platforms and has a limited level of helpful online support. 

The Currently Implemented Processing Scheme 

A simple overview of the current processing flow is as follows: I) SeaWiFS LAC 
telemetry captured at an NRC site and converted immediately to NASA LO formatted file, II) 
Data sent to the satellite data storage facility at Naval Oceanographic Center (NAVO) at Stennis 
Space Center in Mississippi via file transfer protocol (ftp) push, UJ) Data are pushed via ftp to 
NRL Monterey from NAVO, IV) A periodic local cron job searches for receipt of new files, and 
if found: V) SeaWiFS processing script commences to calibrate, clean, apply atmospheric 
corrections, and generate products, VI) Products are uploaded to a password protected NIPRNET 
website. Latency of near real-time products is generally on the order of 1.5 hours from time of 
ingest at the NRC site. The remainder of this section provides additional details on the steps 
outlined above. 

NRC Stations 

The mid-latitude NRCs included in the current NRL processing are Rota [36.63N, 
6.34E], Bahrain: [26.16N, 50.39E], and Yokosuka [35.28N, 139.67E]. Each station captures 
SeaWiFS passes on or around its local noontime hours (1100—1400 Z for Rota, 0700—1000 Z 
for Bahrain, and 0130—0430 Z for Yokosuka). NRL receives the NASA L0 data from NAVO as 
Gnu-compressed files typically 1 hour after the nominal capture time. File sizes vary according 
to pass coverage, with typical values ranging from 20 to 80 Megabytes. The 55 Gigabyte hard 
drive currently used as a storage bay for raw data, algorithms, and products is sufficient for 
several months of online storage. All data are eventually archived to CDROM or tape storage in 
a hierarchical structure partitioned by station, year, month and day. 

Data Calibration 

A cron job launched every half hour initiates the processing of newly received SeaWiFS 
data. A loop over all files present in the ftp "receive directory" checks for any new and 
unprocessed data. If a new file has arrived, it is uncompressed and copied to the local storage 
drive for subsequent processing to LIB calibrated radiances (expressed in mW/cm 2 -sr-jim). 
Binary executables taken from the SeaDAS software package are used for this purpose, and 
require an updated SeaStar ephemeris for correct earth location (SeaWiFS requirement of 1 -pixel 
geolocation accuracy). These data are available at ftp.samoa.g s fc.nasa.gov . and are updated on 
the local system if the current elements file is found to be more than one day old. Any errors 
incurred during the uncompressing or calibration steps are reported to a log file and the script is 
exited gracefully. 

Sector Definition and Noise Corrections 

Upon completion of LIB processing, an DDL script for customization of the data is 
spawned. Latitude, longitude, sun/sensor geometry, and SeaWiFS channel variables are extracted 
from the HDF data. A table of pre-defined sectors for each station is included and readily 
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updated for new operational areas of interest or one-time case studies. Definition of these sectors 
is based on the specification of latitude/longitude boxes (lower left and upper right corners). 
Owing to the scanning geometry, regions near the edge of the swath suffer spatial resolution 
degradation. As a result, imagery for these pre-defined sectors may include only partial coverage 
at compromised quality. A minimum data check of 20% coverage over the sector in question is 
enforced prior to proceeding with the processing. To ensure that the optimal information content 
for every SeaWiFS pass is provided, a "floater" sector is included which selects a box data away 
from the swath edges (poor resolution) and the horizons (avoiding numerous scan line drop-outs). 
The trade off for high image quality is a random coverage area from pass to pass. 




Fig-2 - Raw SeaWiFS data from NRC station (left) and noise/line corrected equivalent (right) 



Once the sector of interest has been defined (or loops over numerous sectors), an embedded loop 
over SeaWiFS channels is performed for corrections of pixel noise and line dropouts. Noisy 
pixels appear in the imagery as a "salt and pepper effect" owing to the presence of spurious very 
high or very low data values. Line dropouts appear in imagery as black scan lines devoid of data 
or filled with noise. These imperfections are present in the data due to a combination of possible 
hardware limitations/problems and clutter associated with nearby structures. Cleaning of the data 
is accomplished in two ways. Lines identified as dropouts are replaced by previous good scan 
lines (bad lines at the very beginning or end of a data stream are replaced by their nearest good 
neighbor). If multiple adjacent dropouts are present, the previously patched line is used once 
again. This correction loses its cosmetic effectiveness for large data gaps, but provides marked 
improvements in cases of spurious line dropouts. Noisy pixels, assumed to be spurious, are 
identified using a threshold based on the absolute values of pixel derivatives computed along a 
scan line. When flagged, they are replaced by an average of the two adjacent pixels on the same 
scan line. Noise pixels on the edges of lines are simply replaced by the value of their nearest 
neighbor. Figure 2 illustrates the result of corrections applied to a subsector of SeaWiFS data 
collected over the Alps. The overall effect of these bad line and pixel noise corrections is 
dramatic, and regarded as an essential step in the production of high-quality SeaWiFS imagery. 
Additional noise corrections (not yet implemented) may be required to correct pixels where one 
of the channels is missing data (but was not flagged by the noise threshold). This results in a 
magenta, yellow, or cyan-biased pixel values in the RGB products. 

Atmospheric Corrections 

The effects of atmospheric molecular (sometimes referred to as "Rayleigh") scattering 
must be removed from the data before creation of true color products. The amount of Rayleigh 
scatter is proportional to X A (wavelength) and hence is strongest for the shorter (e.g., blue, violet) 
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wavelengths. If left uncorrected it will manifest in the imagery as a semi-transparent milky blue 
haze increasing in opacity toward the edges of the swath (as higher sensor angles peer through an 
optically thicker atmospheric path), detracting from image contrast and clarity. To remove this 
undesirable component of the satellite signal, we use a radiative transfer model (Miller et ah, 
[2000]) to simulate the clear sky Rayleigh scatter over a dark surface as a function of the solar 
zenith, satellite zenith, and solar-satellite relative azimuth angles. Pre-calculated look-up tables 
were created for each SeaWiFS channel and the correction was applied on a pixel-by-pixel basis, 
provided the sun/satellite geometry in the Level IB data. The primary caveat associated with this 
correction is the assumption of a cloud-free Rayleigh contribution that in some instances may 
result in an over-correction of pixels containing high clouds near the edges of the swath. Visual 
inspection of imagery has not indicated a significant distortion owing to this effect, however. 
Regions near the edge of the SeaWiFS swath should in general be interpreted with cautious 
regard to the possible distortions owing to larger atmospheric path lengths, three-dimensional 
effects (parallax), and possible bow-tie effects. 

Scaling 

Once the data have been calibrated, cleaned of noise, and corrected for Rayleigh scatter, 
the spectral radiances (I£ are converted to equivalent isotropic reflectances (Rx) according to the 
simple relationship: 

where ^ is the cosine of the solar zenith angle and F 0 is the solar spectral flux at wavelength A.. 
By definition, the isotropic reflectance represents the reflectance produced by an equivalent 
lambertian surface, and varies from 0.0 to 1.0. Any spurious values greater than 1.0 or less than 
0.0 are thresholded by these upper and lower bounds. The reflectances are scaled by log w in an 
attempt to equalize the imagery and prevent the very bright cloud/snow pixels from dominating 
the information content of the products. These scaled reflectances are byte-scaled over the 256- 
element color table range of the IDL palette, with a minimum logarithmic reflectance of -1.65 
mapped to the minimum saturation value (0) and 0.0 mapped to the maximum (255). 

The scaled data are warped to a mercator projection corresponding to the sector of 
interest. This procedure is currendy the most computationally time-consuming component of the 
processing, and options for optimization or alternative approaches are currently being explored. 
A segmentation fault error discovered in the IDL "map patch" command was discovered during 
this development and is currently being investigated by Research Systems Incorporated (RSI). 
Resolving this possible code problem may provide additional computational efficiency to the 
current code. 

The final image products are created as jpegs at 80% of full quality. File sizes range 
from 200 to 800 Kilobytes depending on the sector and image information content. The true 
color, vegetation enhancement, and dust enhancement products for various stations and 
subsectors are copied to a password-protected online holding directory on the NRL Satellite 
Meteorology web page (http://kauai.nrlmrv.navv.mil/sat products.html ). 
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Product Examples 

7>u£ Co/or Enhancement 

To produce the true color images, the blue (412 nm), green (555 nm), and red (670 nm) 
SeaWiFS data are loaded into the respective blue, green, and red IDL "color guns." The data 
comprise a single matrix of dimension (n,m,3) where n=along-scan direction and m=along-track 
direction. The resultant imagery is similar, but not identical, to what would be observed by the 
human eye from the vantage point of outer space; clouds as shades of gray to white, ocean 
surfaces as shades of blue, and land features as earth-tones and green in regions of heavier 
vegetation. 

Examples of true color SeaWiFS imagery are shown for the Rota, Bahrain, and Yokosuka 
regional centers in Figures 3-5, respectively. Figure 3 demonstrates the immediate ease with 
which significant airborne dust is identified over oceans. The image depicts a widespread 
easterly plume of Saharan dust captured within the synoptic flow of a baroclinic system off the 
coast of Spain and Portugal. SeaWiFS is useful in monitoring these dust plumes as they migrate 
across the Atlantic Ocean (often reaching the eastern United States). The green colors in littoral 
regions of the Persian Gulf in Figure 4 correspond to a higher amount of bioactivity (e.g., 
plankton) in nutrient-rich waters. Cloud streets over land and evidence of biomass burning in the 
form of diffuse smoke plumes are apparent over Southeastern Asia in Figure 5. The NRC in 
Yokosuka also captures seasonal dust storms blowing off the Gobi desert. Owing to the general 
circulation of the atmosphere, large Gobi and Saharan dust storms may reach the continental 
United States (the former crossing the Pacific with the "Roaring Fourties", and the latter crossing 
the Atlantic with the Tropical Trades). Implications of these teleconnections are discussed by 
Westphal et of. [1988], Prospero et a/.[1987], and Prospero [1999]. 
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Fig-3 - SeaWiFS HRPT true color image of the Eastern Atlantic Ocean, Iberian Peninsula, and Northwestern Africa, 
captured by NRC in Rota, Spain (02/1 3/01, 1255Z). A large dust plume wrapping into a baroclinic system is indicated 
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SeaWiFS BHR 2001/11/13 08412 Naval Research Lab Mry True Color 
(Red-670nm. Green -555nm. Bue-412nm1 




Fig-4 SeaWiFS HRPT true color image of the Persian Gulf area, captured by NRC in Bahrain. Airborne dust is 
indicated. Phytoplankton is responsible for the green pigmentation in the littoral regions 
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Fig-5 SeaWiFS HRPT true color image of the Korean Peninsula, captured by NRC in Yokosuka Japan. Yellow 

arrows indicate areas of biomass smoke/pollution 



Vegetation Enhancement 

Because green vegetation possesses a relatively high albedo compared to soil in the 
"reflective" or shortwave infrared (0.7-1.3 fim), a combined vis/short wave-IR product that 
enhances vegetation regions is generated in a way similar to the true color enhancement described 
above. This vegetation enhancement product is similar to what has been produced for many years 
for land-usage studies with Landsat TM data. The 670 nm (red) channel in the previously 
discussed RGB true color enhancement is replaced by the 865 nm (SeaWiFS channel 8) channel 
data. This results in green vegetation appearing as bright red and appears usually in sharp 
contrast to surrounding regions devoid of vegetation. An example of this product is shown in 
Figure 6. The vegetation enhancement may prove useful in identifying potential fetch areas for 
airborne dust. 
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Fig-6 - True color (top) and vegetation enhancement (bottom). Bright red indicates regions with abundant plant life. 
Note seasonal changes in Fall vegetation compared against Summer (Figure 5) 
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Airborne Dust Enhancement 

The detection and tracking of airborne dust from satellite has been a long-standing 
interest of the Navy and academia alike, and there is no shortage of papers treating the subject 
(e.g., Shenk and Quran [1974], Carlson [1978], Legrand et al. [1987], Lee [1989], Jankowiak and 
Tanre [1991], Tanre and Legrand [1991]). These earlier works made use of contrast and spectral 
signatures to identify and retrieve properties of dust against land and ocean backgrounds. An 
important capability of SeaWiFS, resulting from its multiple narrowband visible channels and 
indicated in the examples above, is discrimination between airborne dust and clouds/smoke via 
true color methods. Fine-grain particles lifted from sparsely vegetated and dry desert surfaces by 
strong frictional winds (typically greater than 10 m/s, usually by meso- and synoptic-scale 
circulation features) give rise to significant plumes of atmospheric dust. These plumes follow the 
low to mid-level atmospheric circulation and may also become incorporated within synoptic scale 
baroclinic systems (giving rise on occasion to unusually "mud rain" storms that have been 
reported to deposit a coating upon everything in their path). 




Fig-7 - Cartoon illustration of preferential absorption of blue light by dust plumes, resulting in yellow reflected light. 
Clouds scatter more uniformly over the visible wavelengths, resulting in white light 

Perhaps more importantly to Navy interests, airborne dust has significant detrimental 
impacts on slant-range visibilities (with visible optical depths that can exceed 3.0 and 
accompanied by visibilities less than 1 nautical mile) and also poses a serious hazard in terms of 
its potential to damage turbine engines. This is of particular relevance to aircraft carrier 
operations in regions of the world prone to widespread dust outbreaks. Because airborne dust 
retains the earth tone hues (usually light brown or yellow, depending on the source region) of the 
underlying surface, it is readily identifiable in most cases from a true color image, especially for 
optically thicker dust plumes. For regions of less dust (where its effects on slant range visibility 
may still be appreciable), an additional enhancement beyond true color is required. NRL 
Monterey has developed a dust enhancement algorithm for SeaWiFS that attempts to provide this 
capability. Enhancement of dust is predicated on basic spectral properties that characterize its 
appearance at visible wavelengths. As indicated in Figure 7, yellow light (characteristic of many 
desert dust plumes) results from the removal of blue light (via preferential absorption by the dust, 
e.g., Patterson et al. 1977). This spectral contrast is even more dramatic between the blue and 
shortwave infrared channels. This preferential absorption does not occur to the same magnitude 
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in liquid/ice clouds (explaining why clouds appear as shades of gray to white, indicating nearly 
equal reflected intensities of red/green/blue light, in a true color composite). By defining a dust 
enhancement parameter "A": 

A = Q 86S ~ Q 4I2 

where a x are reflectances at wavelength X (nm), any dust component of the image will be 
enhanced. Cloud pixel values, which are relatively spectrally flat, will be diminished in 
brightness owing to the small numerator difference. The logarithm of the A parameter is scaled 
between -0.45 and 0.20 and loaded into the red channel of the RGB composite. 

Owing to the significantly reduced red component for cloud pixels in the dust 
enhancement (small values of the A parameter) image, clouds typically appear as shades of light 
blue or cyan. Retaining the original gray/white cloud hues in the dust enhancement (for cosmetic 
puiposes only) requires the implementation of a simple cloud mask— pixels flagged as cloudy can 
be re-assigned to their original red channel values. Clouds in the imagery are identified first by 
computing the mean scaled magnitude and standard deviation of the 412, 555, and 670 nm 
channels. Mean magnitudes exceeding 50% and standard deviation less than 2.5% are used as 
thresholds for positive cloud identification. The simple premise behind choosing these thresholds 
is that clouds are in general bright (high brightness magnitude) and relatively spectrally flat over 
the visible (low standard deviation). Pixels flagged as "cloudy" are reset to the red channel value 
in the A parameter. As the dust enhancement enhances land surfaces also, a land mask is applied 
to this product. 

^ Figure 8 demonstrates the dust enhancement product covering the Persian Gulf region. 
The top panel pair reveals that while the true color imagery is useful in identifying the significant 
dust plumes, the dust enhancement reveals a far more extensive coverage of airborne dust 
circulating within the Gulf. The center and lower panel pairs further illustrate this enhancement 
capability. While not necessarily posing the same hazard to aircraft engines as the optically thick 
(bright) plumes, these areas of enhanced aerosol have detrimental impacts on slant-range 
visibility. Although more difficult to detect than over the oceans, the presence of significant dust 
over bright and/or sparsely vegetated land surfaces can in some cases be inferred by noting 
contrast reduction in the background terrain. The lack of an infrared channel on the SeaWiFS 
instrument precludes the technique of detecting elevated (cooler) dust above a warm surface 
background via emission measurements. The utility of EOS-MODIS, which provides 36 
channels (including narrowband red, green blue, NIR and LWIR bands), will be explored in this 
capacity. 
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— METOC Post Deployment Report 
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Fig-8 Examples of airborne dust enhancements over the Persian Gulf. Optically thin dust plumes appear as darker , 
tones of red. Land regions in the dust enhancement have been masked 1 
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Other Utilities 

The global coverage afforded by the polar orbiting platform provides SeaWiFS with 
opportunities to survey many atmospheric and terrestrial phenomena. Its utility as an additional 
sensor for tropical storm monitoring (e.g., Typhoon Pabuk shown in Figure 9) may increase in 
relevance in light of the current mechanical problems aboard the Geostationary Meteorological 
Satellite (currently providing coverage of the Tropical Western Pacific and the primary satellite 
observation platform supporting the Joint Typhoon Warning Center). Unexpected events provide 
unique opportunities to examine additional capabilities of the SeaWiFS channels. The eruption of 
Mt. Etna, shown in Figure 10, demonstrates an additional utility of the dust enhancement 
product— revealing features of the ash plume difficult or impossible to detect in panchromatic or 
even true color renditions. Depending on their composition and concentration (not all volcanic 
plumes are comprised of materials that will be enhanced by the currently detailed dust 
enhancement), such plumes may be tracked around the globe, providing invaluable guidance to 
the commercial airline industry for avoiding these potentially serious flight hazards. The 
detection capabilities of SeaWiFS offer a valuable asset to Volcanic Ash Advisory Centers 
(VAACs; e.g., http://www. ssd.noaa.govA^AAC/washington.htmn distributed world-wide whose 
charter is to monitor and track these plumes globally, and also the Geneva-based World 
Meteorological Organization (WMO; e.g., http://www.wmo.ch/indexfl^sh.htmlV which 
coordinates global scientific activity among an international panel of participants. 

Ongoing Developments 

New Processing Formats 

The SeaWiFS processing as outlined in this report is well suited fot the creation and 
analysis of high quality SeaWiFS imagery products in an automated, near real-time framework. 
A Cooperative Research and Development Agreement (CRADA) exists between NRL Monterey 
and the SeaSpace Corporation (whereby a TeraScan license is provided in return for NRL- 
developed algorithms suitable for inclusion to the TeraScan processing suite). Furthermore, the 
NRCs mentioned in this report all use the TeraScan system as their ingest and processing 
platform. 

For the reasons cited above, efforts to reproduce this processing in a fashion entirely 
independent of third party software have resulted in an analogous package that assumes receipt of 
SeaWiFS data in TeraScan Data Format and proceeds to customize the data in the FORTRAN-90 
programming language. True color and dust enhancement image rendering is then completed 
using the "rgbimage" TeraScan intrinsic function. This FORTRAN-based version is actually 
more efficient and hence runs faster than its IDL counterpart. It is currently being tested for 
robustness in anticipation of transition to operational use. Until an HDF-to-TDF converter is 
written that is capable of handling the multiple sensor zenith angle settings characteristic of the 
SeaWiFS sensor (or the data stream changes from L0 to TDF), the default NRL Monterey in- 
house processing will likely remain IDL-based. 

Blended Products 

SeaWiFS is only one of several telemetries available for processing at NRL Monterey. It 
is a straightforward procedure to incorporate data from these other sensors within the same image 
to produce a multi-sensor product. For example, the masked-out land areas in the dust 
enhancement product may be filled in by temporally matched geostationary data (e.g., Meteosat 5 
over Bahrain, MeteoSat 7 over Rota, and GMS over Yokosuka). Infrared information can be 
used to enhance regions of elevated dust over the land in certain situations. As NRL also has at 
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its disposal analysis and forecast data from its mesoscale (Coupled Ocean-Atmosphere Mesoscale 
Prediction System)™ and global scale (Navy Operation Global Atmospheric Prediction System) 
models, it may also prove useful to blend these data into the products. Using again the example 
of the dust enhancement product, inclusion of the COAMPS™ 925 mb wind field in the form of 
vectors overlaid upon the imagery may serve not only as a qualitative validation of the model 
analysis but also provide the METOC officer in the field with a means to predicting the migratory 
behavior of significant dust plumes in the short term. Overlay of short-term forecasts of the 
frictional wind over land (winds near the surface that are responsible for lifting the dust into the 
atmosphere) may also be of use. Of course, any such multi-sensor and model-fusion products 
will require a training module identifying the various components, strengths, and weaknesses. 

These data are of immediate relevance to both the initialization and validation of the 
Navy Aerosol Analysis and Prediction System (NAAPS; a global, multi-component analysis and 
modeling capability to produce aerosol products from satellite data and surface-based 
measurements). The goal of NAAPS is to improve the forecasting/nowcasting of visibility 
restrictions caused by aerosols, including dust, smoke, and volcanic ash. ' Information and 
forecast products associated with NAAPS may be found at http://www.nrlmrv.navv.mil/aemsol 

The image products discussed in this report are intended for use in the field by Navy 
METOC officers for at-a-glance identification of dust storms and/or other features of interest; 
they are qualitative in nature. In other words, these products should not be misinterpreted as 
providing a quantitative flag for delineating the scene for dust or other constituents. Quantitative 
aerosol optical properties and associated slant-range visibility products are retrievable from 
SeaWiFS channels. Although the topic of aerosol retrievals is beyond the scope of the current 
work, extensive work in this area by several colleagues at NRL Monterey (D. Westphal, P. 
Flatau, K. Richardson, A. Kuciauskas) and the Naval Postgraduate School (P. Durkee) is 
ongoing. 

Extension to MODIS Processing 

The MODIS instrument offers a similar true color imaging capability together with 
improved spatial resolution (up to a factor of 4 over SeaWiFS HRPT). Figure 11 provides an 
example of the MODIS dust enhancement over the Northern Arabian Sea. Plumes of dust over 
ocean and dry lakebeds on the desert interior (often serving as sources for dust) are enhanced as 
shades of pink. Of specific relevance to the dust enhancement product is the availability of 
thermal infrared channels that enable the discrimination of elevated (cooler) dust plumes over 
warm backgrounds. On-board calibration of MODIS radiances supports accurate physical 
retrievals of aerosol optical properties. A technique for producing 250m true color imagery using 
500m blue/green and 250m red channels has also been developed. This high spatial resolution 
capability is enhanced by temporal coverage when considering the Terra (MODIS: 1030 
descending node), SeaStar (SeaWiFS: 1200 descending node), and Aqua (MODIS: 1330 
ascending node; platform to be launched in April 2002) platforms in concert. This quasi- 
geostationary 3-hour loop will provide valuable information on cloud and/or dust motion. 
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WHILE OPERATING IN THE MAS IN SUPPORT OF OPERATION ENCORING FREEDOM ( OEF ) , THE PRIMARY METOC- 
RELATED IMPACT TO OPERATIONS WAS DECREASED VISIBILITY IN N, W, AND S AFGHANISTAN DOE TO STRONG 
WINDS U] IN ONE CASE, THE EXTENT OF THE SUSPENDED DOST AREA RANGED WELL OUT INTO THE NORTHERN 
ARABIAN SEA, WITH VISIBILITY LESS THAN 1 NM. 
— METOC Post Deployment Report 



Fig- 11 - MODIS dust enhancement (0.5 km resolution) over the Northern Arabian Sea, true color over land. Dust 
over the ocean and dry lake beds appear as shades of pink. Accompanying text from Navy MetOc post deployment 
report corresponding to the storm observed in this image 



Additional sensors with similar true color imaging capabilities (e.g., GLI and POLDER aboard 
ADEOS n, MISR aboard Terra, OSMI aboard KOMPSAT, OCM aboard OCEANS AT, FY-1C, 
TM/ETM aboard the LandSat series, etc.) would in some cases provide additional temporal 
information. The potential METOC benefit to the war fighter in terms of from these state of the 
art instruments cannot be understated. 
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Conclusion 

A new processing capability for generation of true color, vegetation enhancement, and 
dust enhancement imagery from the 8-channel SeaWiFS instrument has been developed by the 
NRL Monterey Satellite Meteorology Applications Section team. Near real-time processing of 
telemetry captured from Navy Regional Centers in Rota, Bahrain, and Yokosuka currently 
populate the following password-protected website on the NRL Monterey satellite page at 1-2 
hour latency: 

http://kauai . nrlmrv.navv.mi 1/archdat/swf nrc/ 

A number of geo-fixed subsectors spanning the coverage zones for each of the centers have been 
defined, in addition to a "floater" sector that roams with the nadir-most pixels of any given 
SeaWiFS pass. Owing to the sensitivity of regions covered by these regional centers and the 
existing proprietary agreement with Orbimage, access to these products has been restricted to 
Navy METOC operational purposes. A preliminary agreement has been reached with Fleet 
Numerical Meteorology and Oceanography Center (FNMOC) to host the products upon their 
secure Internet website (SIPRNET) while NRL awaits procurement of its own SIPRNET 
capability. A number of other hosting options that would provide 24/7 support of the SeaWiFS 
product suite via secure protocol are currently being explored. 

The relevance of high-resolution imagery to Navy operations in the context of 
aerosol/dust detection is evidenced by numerous METOC post-deployment reports. An excerpt 
from one such report, speaking to METOC-related problems affecting deployments during 
Operation Enduring Freedom (Afghanistan), is included below: 

"On several occasions weapons laser lock was difficult due to blowing dust and 
sand, or low ceilings. Occurred primarily over N, C, and W Afghanistan. [...] successfully 
used high-res satellite imagery for short-term forecasts of blowing dust events in the AG. 
The three blowing dust events that resulted in mission cancellations were all accurately 
forecast" 

-LCDR K. Schwingshakl, USS Enterprise, METOC Post Deployment Report 

Perhaps the greatest testimony to the potential utility of these products has been provided in the 
form of feedback from METOC officers who have seen these products first-hand. The following 
recently received comments from METOC officers with specific regard to the NRL SeaWiFS 
products outlined in this report indicate a strong interest in obtaining real-time imagery onboard 
earners for use in both plenary and tactical guidance: 



"Imagery is fantastic — and very operationally useful" 
-LCDR Ray Chartier, NSAWC METOC Officer, NPMOD Fallon 

"Are these images on SBPR??? If so where...great stuff." 
-LCDR Mike Angove, Staff Meteorologjst/Oceanographer, Carrier Group Three 

Two main points become clear from the above commentary: 1) SeaWiFS and other sensors 
offering true color imaging capability are of significant value and in high demand for Navy 
operations, and 2) in spite of the fact that such data has been available for some time now from 
the satellite platform (e.g., LandSat-TM) this capability has not been transitioned effectively to 
the end-users who require it (i.e., those deployed in the field). In turn, the twofold intent of this 
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report has been to demonstrate the current in-house processing capabilities for the SeaWiFS 
sensor and articulate this operational demand. A logical extension to these statements is that 
there exists a pressing need for X-band receiving stations in perennially sensitive regions of the 
world for the purpose of capturing MODIS and other high-data-rate (both multi- and hyper- 
spectral) satellite telemetries. The potential savings in terms of both loss of life and equipment 
(e.g., Cantu [2001] reports a 51 million dollars annual loss in Navy equipment attributed to 
visibility-related meteorological phenomena such as aerosols, dust, and atmospheric trapping 
layers) owing to improved METOC guidance gained from these state-of-the-art sensors far 
outweighs the one-time cost of the receiving stations themselves. 
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# Call external software to process data (clean, remove rayleigh, scale) 
print "Calling external processing sof tware. \n"; 

system ( n process_modis n ) == 0 or die "Problem with process_modis : $!"; 

##SNIPPET FROM process_modis_color . f 90 PERTAINING TO DE** 
#! 1) Atmospheric Correction (snippet) 

# airmass = 1.0/mu_sun + 1.0/mu_sat 

# trans_sat = exp (-1 . 0*ray_taus (i) /mu_sat) 

# trans_sun = exp (-1 .0*ray_taus (i) /mu_sun) 

# tro - proc_chandat*pi/ ( F_sun_adj (i) * mu_sun * exp (-1 . O*airmass*uoz*aoz (i) ) 
) 

# trayu = ( (2.0/3.0 + mu_sat ) + (2 . 0/3 . 0 - mu_sat ) *trans_sat ) / & 

# ( 4.0/3.0 + ray_taus(i) )• 

# trayd = ( (2.0/3.0 + mu_sun)+(2.0/3.0 - mu_sun) *trans_sun) 7 & 

# ( 4.0/3.0 + ray__taus(i) ) 

# tro = (tro-ray_ref ) / ( trayu*trayd*exp (-1 .0*airmass* (ao2 (i) + awv(i))) ) 
#' proc_chandat = tro/ ( 1 . 0+raysa (i) *tro) ! units of reflectance [0-1] 

#! 2) Truncate negative numbers to 0, if they occur 

#! Truncate numbers greater than unity to unity, if they occur 

# WHERE (proc_chandat < 0.0223) proc_chandat = 0.0223 

# WHERE (proc_chandat > 1.0) proc_chandat =1.0 
#! 3) Scale the processed data if applicable 

# print *, "Scaling the data by log_10" 

# proc_chandat = LOG10 (proc_chandat) 
# ! 4 ) Dust Enhancement 

# print *, "Creating dust enhancement" 

# WHERE ( ( (chan03_dat+chan02_dat) > 0.0) .AND. (chan02_dat > 0.0223) ) 

# dust_enhance = (chan02_dat-chan03_dat ) / (chan02_dat+chan03_dat ) 

# ELSEWHERE 

# dust_enhance = -0.40 

# ENDWHERE 

#! 5) Truncate between desired values 

# WHERE (dust_enhance < -0.40) dust_enhance = -0.40 

# WHERE (dust_enhance > 0.15) dust_enhance =0.15 
###END SNIPPET### 

# 

# Dust Enhancement 

# 

if ($have_ir) { 

print "Creating a \"dust enhanced\" image: ${OUTNAME} .dustlKM. $SECTOR.x. tif f \n"; 
### 

# Water component 
### 

print "Doing ocean component ... \n"; 

$min_blue_scl = -1.45; 
$max_blue_scl =0.0; 
$min_green_scl = -1.45; 
$max_green_scl = 0.0; 
$min_de_scl = -0.40; 
$max_de_scl =0.15; 

system ("emathp exprjvars= f ref_4 69 ref_555 DE f num_exprs=3 
"yl_expr= f land ? $min_blue_scl : xl' ". 
"y2_expr= , land ? $min_green_scl : x2 f ". 
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"y3_expr= , land ? $min_de_scl : x3' 
"use_degrees= poly_size=100 box_sides='l 1\ 

"min_good= ' 1 1 save_exprs= f 1 2 3' var_names='Blue_ocean Green_ocean 
Red_ocean 1 " . 

"var_units= var_types=' float float float 1 scale__of f sets= scale_f actors= 
"reflec_assem_$pid tmp_ocean_$pid") == 0 or die "Problem with emathp: 

$!"; 

system ("imscale include_vars= , Blue_ocean' image_colors=250 bg__color=250 
bad_color=0 

"over_color=249 under_color=0 max_width=$SAMPLES max_height=$LINES 
fixed_£ize=yes 

"zoom_factor=l real_resample=no est_range=no min_value=$min_blue_scl 

max_value=$max_blue_scl invert_scale=no ■ 

,f north_up=yes tmp_ocean_$pid tmp_blue_$pid") == 0 or die 
"Problem occurred during imscale of blue ocean-component of dust 

enhancement : $ ! " ; 

system ("expim image_f ormat=tif f image_var=Blue_ocean image_colors=250 
color_palette=black-white " . 

"draw_indexes= draw_names= tmp__blue_$pid blue_ocean.tif f_$pid") == 0 or 

die 

"Problem occurred during expim of blue_ocean. tif f_$pid creation: $!"; 

system ("imscale include_vars= 1 Green_ocean ' image_colors=250 bg_color=250 
bad_color=0 ". 

"over_color=249 under_color=0 max_width=$ SAMPLES max_height=$LINES 
fixed_size=yes ". 

"zoom_factor=l real_resample=no est_range=no min_value=$min_green_scl 

max_value=$max_green_scl invert_scale=no 

"north_up=yes tmp_ocean_$pid tmp_green_$pid" ) == 0 or die 
"Problem occurred during imscale of green ocean-component of dust 

enhancement : $ ! " ; 

system ("expim image_f ormat=tif f image_var=Green_ocean image__colors=250 
color_palette=black-white " . 

"draw_indexes= draw_names= tmp_green_$pid green_ocean. tiff_$pid") == 0 

or die 

"Problem occurred during expim of green_ocean.tif f_$pid creation: $!"; 

system ("imscale include_vars=' Red_ocean' image_colors=250 bg_color=250 
bad_color=0 ". 

"over_color=24 9 under_color=0 max_width=$ SAMPLES max_height=$LINES 
f ixed_size=yes ". 

"zoom_factor=l real_resample=no est_range=no min_value=$min_de_scl 
max_value=$max_de_scl invert_scale=no ". 

"north_up=yes tmp_ocean_$pid tmp_de_$pid") == 0 or die 

"Problem occurred during imscale of red (DE) ocean-component of dust 
enhancement : $ ! " ; 

system ("expim image_format=tif f image_var=Red_ocean image_colors=250 
color_palette=black-white " . 

"draw_indexes= draw_names= tmp_de_$pid red_ocean*tiff_$pid") == 0 or die 
"Problem occurred during expim of red_ocean.tif f_$pid creation: $!"; 

# Make RGB image 

system ("tifftopnm red_ocean. tif f_$pid > red_ocean.ppm_$pid") == 0 or die 
"Problem with tifftopnm call: $!"; 
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system ("tifftopnm green_ocean . tif f_$pid > green_ocean .ppm_$pid" ) == 0 or die 
"Problem with tifftopnm call: $!"; 

system ("tifftopnm blue_ocean.tif f_$pid > blue_ocean.ppm_$pid") == 0 or die 
"Problem with tifftopnm call: $!"; 

system ("ppmtopgm red_ocean.ppm_$pid > red_ocean.pgm_$pid") =?= 0 or die 
"Problem with ppmtopgm call: $!"; 

system ("ppmtopgm green_ocean.ppm_$pid > green_ocean.pgm_$pid" ) — 0 or die 
"Problem with ppmtopgm call: $!"; 

system ("ppmtopgm blue_ocean.ppm_$pid > blue_ocean.pgm_$pid") == 0 or die 
"Problem with ppmtopgm call: $!"; 

system ("rgb3toppm red_ocean.pgm_$pid green_ocean .pgm_$pid blue_ocean.pgm_$pid 
> tmp_ocean.ppm_$pid") == 0 or die "Problem with rgb3toppm call: $!"; 

system ("convert ppm: tmp_ocean.ppm_$pid tiff : tmp_ocean_one . tif f__$pid") == 0 
or die "Problem occurred during convert of modis dust enhancement 
product : $ ! " ; 1 

### 

# Sunglint mask 
### 

print "Making sunglint mask for use with dust enhancement imagery ... \n"; 
$min_sunglint = 30.0; 
$delta = 0.25; 

system ("emathp expr_vars=' sunglint ref_853' num_exprs=l ". 

"yl_expr=' land ? badval : 
(xl<($min_sunglint+$delta) ) && (xl> ($min_sunglint-$delta) ) ? ". 

"100 : (xl<($min_sunglint-$delta) ) ? x2 : badval 1 ". 
"use_degrees= poly_size=100 box_sides =, l 1' ". 
"min_good= f 1 1 save__exprs= 1 1 1 var_names= 1 glintmask ' " . 
"var_units= var^ypes^ 1 float ' scale__of f sets= scale_factors= ". 
"reflec_assem_$pid tmp_glintl_$pid") == 0 or die "Problem with emathp: 

$!"; 

system ("imscale include_vars= f glintmask' image_colors=250 bg_color=250 
bad_color=0 ". 

"over_color=252 under_color=2 max_width=$ SAMPLES max_height=$LINES 
f ixed_size=yes 

"zoom_factor=l real_resample=no est__range=no min__value=$min_green_scl 

max_value=$max_green_scl invert_scale=no " . 

"north_up=yes tmp_glintl_$pid tmp_glint_$pid") == 0 or die 
"Problem occurred during imscale of green ocean-component of dust 

enhancement : $ ! " ; 

system ("expim image_format=tif f image_var=glintmask image_colors=250 
color_jpalette=black-white " . 

"draw_indexes= ' 252 1 draw_names= 1 yellow 1 tmpjglint_$pid 
tmp_glint.tiff_$pid") == 0 or die 

"Problem occurred during expim of tmp__glint . tiff_$pid creation: $!"; 

# Merge Sunglint Mask and Ocean image 

system ( "imgmerge background_f ile=tmp_ocean_one . t if f_$pid 
foreground_file=tmp_glint . tiff_$pid clear_color=black " . 

"output_format=tif f tmp_ocean.tif f_$pid") — 0 or die "Problem with 
imgmerge: $!"; 

### 

# Land component 
### 
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print "Doing land component- . . \n"; 
#.### 

$min_blue_scl = -1.45; 
$max_blue_scl = 0.0; 
$min__green_scl = -1.45; 
$max_green_scl = 0.0; 
$min_red_scl = 1.30; 
$max_red_scl = 2.70; 
$min_rmblnd = -0.10; 
$max_rmblnd = 0.25; 

@stat = 'stats include_vars=modis_ch31 printout=no ref lec_assem_$pid x ; 
@stat2 - split (/ +/,$stat[4] ) ; 
$min_31 = $stat2[3]; 
$max_31 = $stat2[4]; 

print "Modis ch31 has min value of $min_31 and max of $max_31\n"; 
$max_irlnd = $max_31; 
if ($max_irlnd < 28.0) { 
$min_irlnd = $max_irlnd-21 . 0; 
} else { 

$min_irlnd = $max_irlnd/4 . 0; 

} 

print "...setting min_irlnd = $min_irlnd, max_irlnd = $max_irlnd\n" ; 
# Should make fmf limits dynamic 

$min_fmflnd = -2.0; 

$max_fmflnd = 2.0; 

$mi n_rmb_i r_di f f = -1.5; 

$max_rmb_ir_dif f = 0.25; 

$min_ch26 = 5.0; 
#### 

system ("emathp expr_vars=' ref_469 ref_555 ref_645 ref_853 modis_ch26 
modis_ch31 modis_ch32 f num_exprs=10 ". 

,f yl_expr='land && ! bad(x6)? pow(10 / x3j -pow(10, xl) :badval' ". 

"y2_expr='land && ! bad(x6)? pow(10,x2) -pow(10,x3) ibadval 1 ". 

"y3_expr= f land && ! bad(x6)? yl-y2 rbadval' ". 

"y4_expr='land && ! bad(x6)? (x6>$max_irlnd) ? $max_irlnd : 
(x6<$min_irlnd) ? $min_irlnd : x6 : badval' ". 

"y5_expr='land && ! bad(x6)? (y4-$min_irlnd) / ($max_irlnd-$min_irlnd) : 
badval 1 ". 

"y6_expr=»land && ! bad(x6)? ((y3-y5)- 
$min_rmb_ir_dif f ) / ($max_rmb_ir_dif f-$min_rmb_ir_dif f ) : badval 1 " . 

"yT^expr^'xV-xG 1 ". 

"y8_expr='land && ! bad(x6)? (y7>$max_fmf lnd) ? $max_fmflnd : 
(y7<$min_fmflnd) ? $min_fmflnd : y7 : badval 1 ". 

"y9_expr= ' (y8-$min_fmf lnd) / ($max_fmf lnd-$min_fmf lnd) 1 " . 

"yl0_expr='land && ! bad(x6)? (x5>$min_ch26) ? 1.0 : 0.0 : badval 1 ". 

"use_degrees= poly_size=100 box_sides='l 1' ". 

"min_good= , l l save_exprs=' 1 23456789 10 ' ". 

"var_names=' red_blu_dif f grn_red_diff yl_y2_diff ir_trunc ir_norm 
y3_y5_diff fmf fmf_trunc fmf_norm wv26_norm' ". 

"var_units= var_types=' float float float float float float float float 
float float 1 scale_of fsets= scale_f actors= " . 

"ref lec_assem_$pid rmb_ir_landdif f_$pid") == 0 or die "Problem with 
emathp: $!"; 

system ("copyvar include_vars= , y3_y5_dif f fmf_norm wv26_norm ir_norm' 
overwrite vars=n " . 
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"rmb_ir_landdiff_$pid ref lec_assem_$pid") == 0 or die "Problem with 
copyvar: $!"; 

$min_topo = 3000.0; 
$topo_scale = 3000.0; 

system ("emathp expr_vars='ref__4 69 ref_555 y3_y5_diff fmf_norm wv26_norm 
modis_ch03 ir_norm elevation' num_exprs=3 " . 

"yl_expr='land && ! bad(x6) ? xl : badval' 

"y2_expr='land && ! bad(x6) ? x2 : badval 1 

"y3_expr='land && ! bad(x6) ? (x8>$min_topo) ? x3+x4-x5+ (1 . 0-x7 ) - 
x8/$topo_scale: x3+x4-x5+ (1. 0-x7) : badval' 

"use_degrees= poly_size=100 box_sides='l 1' 

n min_good= 1 1 1 save_exprs=' 1 2 3' var__names=' Blue land Green land 
Red_land' 

"var_units= var_types=' float float float 1 scale offsets= scale factors^ 
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"reflec_assem_$pid tmp_land $pid") == 0 or die "Problem with emathp: 

$!."; 

system ("emathp expr_vars='Red_land' num_exprs=l 
"yl_expr='land ? xl : $min_red_scl ' ". 
"use_degrees= poly_size=100 box_sides='l 1' " . 
"min_good= ' 1 ■ save_exprs= 1 1 ' var__names= 1 RedJLand ' " . 
"var_units= var_types=' float 1 scale_of f sets= scale_factors= ". 
"tmp_land_$pid tmp_land red $pid") == 0 or die "Problem with emathp: 
$!";"" 

system ("copyvar include_vars= f Red_land' overwrite_vars=yes " . 

"tmp_land_red_$pid tmp_land $pid") == 0 or die "Problem with copyvar: 

$!"; 

system ("imscale include_vars='Blue_land' image_colors=250 bg color=250 
bad_color=0 ". 

"over_color=24 9 under_color=0 max_width=$SAMPLES max_height=$LINES 
f ixed_size=yes ". 

"zoom_factor=l real_resample=no est_range=no min_value=$min_blue_scl 

max_value=$maxjDlue_scl invert_scale=no " . 

"north_up=yes tmp_land_$pid tmp_landblue_$pid") == 0 or die 
"Problem occurred during imscale of blue land-component of dust 

enhancement: $!"; 

system ("expim image_f ormat=tif f image_var=Blue_land image_colors=250 
color_palette=black-white ". 

"draw_indexes= draw_names= tmp_landblue_$pid blue_land. tiff_$pid") == 0 

or die 

"Problem occurred during expim of blue_land. tif f_$pid creation: $!"; 

system ("imscale include_vars='Green_land' image_colors=250 bg_color=250 
bad_color= ". 

"over_color=24 9 under_color=0 max_width=$SAMPLES max_height=$LINES 
f ixed_size=yes ". 

"zoom_factor=l real_resample=no est_range=no min_value=$min_green_scl 

max_value=$max_green_scl invert_scale=no ". 

"north_up=yes tmp_land_$pid tmp_landgreen_$pid") == 0 or die 
"Problem occurred during imscale of green land-component of dust 

enhancement : $ ! " ; 

system ("expim image_format=tif f image_var=Green_land image_colors=250 
colorjpalette=black-white 
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"draw_indexes= draw_names= tmp_landgreen_$pid green_land. tif f_$pid") == 

0 or die 

"Problem occurred during expim of green_land. tif f_$pid creation: $!"; 

system ("imscale include__vars= f Red_land' image_colors=250 bg_color=250 
bad_color= 

"over_color=24 9 under_color=0 max_width=$ SAMPLES max_height=$LINES 
fixed_size=yes 

"zoom_factor=l real_resample=no est_range=no min_value=$min_red scl 
max_value=$max_red__scl invert_scale=no ~~ 
"north_up=yes tmp_land_$pid tmp_landred_$pid") == 0 or die 
"Problem occurred during imscale of red land-component of dust 
enhancement : $ ! " ; 

system ("expim image_format=tif f image_var=Red_land image_colors=250 
color_palette=black-white " . 

"draw_indexes= draw_names= tmp_landred_$pid red_land. tif f_$pid") == 0 or 

die 

"Problem occurred during expim of red_land. tif f_$pid creation: $!"; 

### 

# Make RGB image 
### 

system ("tifftopnm red_land. tif f_$pid > red_land.ppm_$pid" ) == 0 or die 
"Problem with tifftopnm call: $!"; 

system ("tifftopnm green_land.tif f_$pid > green_land.ppm_$pid" ) == 0 or die 
"Problem with tifftopnm call: $!"; 

system ("tifftopnm blue_land. tif f_$pid > blue_land.ppm_$pid") == 0 or die 
"Problem with tifftopnm call: $!"; 

system ("ppmtopgm red_land.ppm_$pid > red_land.pgm_$pid") == 0 or die "Problem 
with ppmtopgm call: $!"; 

system ("ppmtopgm green_land.ppm_$pid > green_land.pgm_$pid") == 0 or die 
"Problem with ppmtopgm call: $!"; 

system ("ppmtopgm blue_land.ppm_$pid > blue_land.pgm_$pid") == 0 or die 
"Problem with ppmtopgm call: $!"; 

system ("rgb3toppm red_land.pgm_$pid green_land.pgm_$pid blue_land.pgm_$pid > 
tmp_land.ppm_$pid") == 0 

or die "Problem with rgb3toppm call: $!"; 

system ("convert ppm: tmp_land.ppm_$pid tiff : tmp_land.tif f_$pid") == 0 

or die "Problem occurred during convert of modis dust enhancement 
product: $!"; 

# Merge Ocean and Land RGB imagery, and metafile image 
system ("imgmerge background_f ile=tmp_ocean. tif f_$pid 

foreground_file=tmp_land. tif f_$pid clear_color=black 

"output_format=tiff ${OUTNAME} .dust_tmp.tif f_$pid") == 0 or die "Problem 
with imgmerge: $!"; 

system ("legend line_l=\"$ { satname} -MODIS $mm/$dd/$yyyy ${hhmm}-Z 1km Dust NRL 
Monterey \" ". 

"line_2=\"Pink=Possible Dust, Black/White=$ {min_sunglint } -Degree 
Sunglint Zone\" line_3= text_height=$TEXT_HEIGHT text_type=bold line_of f set=0 
sample_of fset=0 ". 

"full_width=no center_text=no outline=yes solid bg=yes $ MASTER $LEGEND") 

0 

or die "Problem with legend: $!"; 
system ("xinject bg_var_name= f master 1 colors_per=2 meta_colors=' 5 0 3 0 5 4' 



"line_widths='2 1 2 1 line_types= marker_types= marker_sizes= 
"image_colors=4 $ MASTER $GRID $COAST $ LEGEND meta . tdf_$pid" ) ==0 
or die "Problem with xinject: $!"; 

system ("expim image_format=tif f image_var=master image_colors=4 
"color_palette=black-white draw_indexes=\"3 4 5 6\" 

"draw_names=\ "violet white blue salmon\" meta. tdf_$pid meta. tif f_$pid") 

== 0 

or die "Problem with expim: $!"; 

print "Merging composite with meta image... \n"; 

system ("imgmerge background_f ile=$ {OUTNAME} .dust_tmp.tif f_$pid 
foreground_file=meta. tif f_$pid clear_color=black ". 

"output_format=tiff ${OUTNAME} . dustlKM. $ SECTOR. x. tif f_$pid") == 0 or die 
"Problem with imgmerge: $!"; - 

# Adjust image quality 
system ("convert -quality 80 tiff :${OUTNAME} .dustlKM.${ SECTOR) .x. tiff_$pid ". 
"jpeg:$ {OUTNAME} . dust 1KM. ${ SECTOR} .x. jpg_$pid") == 0 " 
or die "Problem occurred during convert of modis dust enhancement 

product : $ ! " ; 



